AFRL-RY-WP-TP-2009-1304 


EFFECT  OF  BANDWIDTH  ON  WIDEBAND-STAP 
PERFORMANCE  (PREPRINT) 

Ke  Yong  Li,  Unnikrishna  S.  Pillai,  and  Joseph  R.  Guerci 
C  &  P  Technologies,  Inc. 


OCTOBER  2007 


Approved  for  public  release;  distribution  unlimited. 

See  additional  restrictions  described  on  inside  pages 


STINFO  COPY 


AIR  FORCE  RESEARCH  LABORATORY 
SENSORS  DIRECTORATE 

WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7320 
AIR  FORCE  MATERIEL  COMMAND 
UNITED  STATES  AIR  FORCE 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite 
1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YY)  2.  REPORT  TYPE 

October  2007  Conference  Paper  Preprint 


4.  TITLE  AND  SUBTITLE 

EFFECT  OF  BANDWIDTH  ON  WIDEBAND- STAP  PERFORMANCE  (PREPRINT) 


6.  AUTHOR(S) 

Ke  Yong  Li  (C  &  P  Technologies,  Inc.) 
Unnikrishna  S.  Pillai  (Polytechnic  University) 
Joseph  R.  Guerci  (Guerci  Consulting) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

C  &  P  Technologies,  Inc.  Polytechnic  University 

317  Harrington  Avenue  Brooklyn,  NY  11201 

Suites  9  &  1 0  """  “  . 

Closter,  NJ  07624-1911 

Arlington,  VA  22207 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


3.  DATES  COVERED  (From  -  To) 

25  May  2006  -  19  October  2007 


5a.  CONTRACT  NUMBER 

FA8750-06-C-01 17 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT 
NUMBER 

62204F _ 

5d.  PROJECT  NUMBER 

5017 _ 

5e.  TASK  NUMBER 

RL 


5f.  WORK  UNIT  NUMBER 

517R151 1 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSORING/MONITORING 
AGENCY  ACRONYM(S) 

AFRL/RYRT 

11.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER(S) 

AFRL-RY  -  WP-TP-2009- 1 3  04 


Air  Force  Research  Laboratory 
Sensors  Directorate 

Wright-Patterson  Air  Force  Base,  OH  45433-7320 
Air  Force  Materiel  Command 
United  States  Air  Force 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 


13.  SUPPLEMENTARY  NOTES 

Conference  paper  to  be  published  in  the  Proceedings  of  the  41st  Annual  Asilomar  Conference  on  Signals,  Systems,  and  Computers, 
held  November  04  -  07,  2007  at  the  Asilomar  Hotel  and  Conference  Grounds,  Pacific  Grove,  CA.  PAO  Case  Number:  WPAFB  07- 
0330;  Clearance  date:  05  Nov  2007.  Paper  contains  color.  See  also,  AFRL-RY-WP-TP-2009-1303  for  a  briefing  chart  version,  and 
AFRL-RY -WP-TP-2009-1305  for  a  postprint  version  of  the  paper. 


14.  ABSTRACT 

A  wideband  signal  occupies  a  finite  bandwidth  that  is  significant  compared  to  its  carrier  frequency.  As  a  result,  when 
transmitted,  its  returns  cause  bandwidth  dispersion  across  the  antenna.  It  is  shown  here  that  the  effect  of  the  finite 
bandwidth  is  to  introduce  a  set  of  uncorrelated  return  signals  for  every  physical  scatter  in  the  field.  Further,  each  such 
uncorrelated  return  contains  a  set  of  coherent  signals  with  different  directional  and  Doppler  components  that  result  from 
a  jittering  effect  both  in  angle  and  Doppler  domain.  As  a  result,  adaptive  clutter  cancellation  using  traditional  processing 
schemes  does  not  work  well.  Although  in  principle  it  is  possible  to  correct  these  decorrelating  effects  by  3D  spacetime 
adaptive  processing  (STAP),  the  present  day  methods  are  quite  costly  and  difficult  to  implement.  In  addition  to  the  new 
wide  band  signal  modeling  framework  mentioned  above,  we  outline  a  hierarchical  processing  scheme  which  has  the 
potential  for  dramatically  reducing  both  processing  and  sample  support  burdens. 

15.  SUBJECT  TERMS 

wideband  space-time  adaptive  processing 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

17.  LIMITATION 
OF  ABSTRACT: 

SAR 


18.  NUMBER 
OF  PAGES 

10 


19a.  NAME  OF  RESPONSIBLE  PERSON  (Monitor) 
Michael  J.  Callahan 

19b.  TELEPHONE  NUMBER  (Include  Area  Code) 
N/A 


1 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


Effect  of  Bandwidth  on  Wideband-STAP 

Performance 


S.  Unnikrishna  Pillai1,  Ke  Yong  Li 1  and  Joseph  R.  Guerci3 

1  Electrical  Engineering  Department,  Polytechnic  University,  Brooklyn,  NY  11201,  pillai@hora.poly.edu 
2C  &  P  Technologies,  Inc.,  Closter,  NJ  07624,  kli@cptni.com 
3Guerci  Consulting,  Arlington,  VA,  irguerci@ieee.org 


Abstract-A  wideband  signal  occupies  a  finite  bandwidth  that 
is  significant  compared  to  its  carrier  frequency  inasmuch  as 
when  transmitted,  its  returns  cause  bandwidth  dispersion 
across  the  antenna.  It  is  shown  here  that  the  effect  of  the  finite 
bandwidth  is  to  introduce  a  set  of  uneorrelated  return  signals 
for  every  physical  scatter  in  the  field.  Further,  each  such 
uneorrelated  return  contains  a  set  of  coherent  signals  with 
different  directional  and  Doppler  components  that  result  from 
a  jittering  effect  both  in  angle  and  Doppler  domain.  As  a 
result  adaptive  clutter  cancellation  using  traditional  processing 
schemes  don’t  work  well.  Though  in  principle  it  is  possible  to 
correct  these  decorrelating  effects  by  3D  space-time  adaptive 
processing  (STAP),  the  present  day  methods  are  quite  costly 
and  difficult  to  implement.  In  addition  to  the  new  wideband 
signal  modeling  framework  mentioned  above,  we  discuss  a  new 
hierarchical  processing  scheme  which  has  the  potential  for 
dramatically  reducing  both  processing  and  sample  support 
burdens. 


b(o)dt)  =  [  1  enm *  •••  (4) 

representing  the  temporal  steering  vector  along  the  Doppler 
frequency 

2VTsm6, 

CO,  = - -  .  PI 

dt  XU 

In  (5),  V  represents  the  platform  speed,  T  the  pulse 
repetition  interval  and  X  the  operating  wavelength. 

In  an  uneorrelated  clutter  and  noise  scene,  from  (1),  we 
get  the  space-time  covariance  matrix  to  be 


I.  Introduction 

In  the  context  of  target  detection  in  a  wideband 
transmit/receive  environment,  consider  an  N  element,  M 
pulse  array  receiving  signals  from  its  field  of  view.  In  the 
narrowband  setup  with  e'"’"'  representing  the  transmit  signal, 
after  down  converting  to  base  band  the  array  output  space- 
time  data  vector  \(t)  from  any  range  bin  has  the  form 

x(0  =  Xa*s^-’ @c)  +  n( 0  (1) 

k 

where  ak  represents  the  random  scatter  return  from  the  kth 
bin  in  the  cross-range  (azimuth)  domain.  Here 

sk  =  s(dk ,  codi )  =  b(codi )  0  a(9k)  (2) 

represents  the  MN  x  1  space-time  steering  vector  with 
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and  its  properties  are  well  documented  [1].  Because  of  the 
angle-Doppler  linear  relationship  in  (5),  in  the  narrowband 
case  for  a  specific  look  direction,  the  adaptive  space-time 
processor  using  (6)  will  null  out  clutter  components  induced 
along  that  direction  thereby  making  target  detection  possible. 

II.  WlDBAND  STAP  PERFORMANCE 

To  obtain  a  similar  understanding  for  the  structure  of  the 
data  and  array  output  covariance  matrix  in  the  wideband 
case,  assume  that  matched  filtering  and  adaptive  processing 
is  employed  in  that  case  as  well.  Consider  a  wideband 
signal  /(f)  with  bandwidth  Ba  that  is  modulated  by 

frequency  coo  to  generate  f(t)ejmJ  which  acts  as  the 

transmit  signal.  When  a  wideband  waveform  is  employed 
in  the  context  of  space  time  adaptive  processing  (STAP), 
the  finite  bandwidth  dispersion  across  the  antenna  must  be 
taken  into  account.  To  understand  the  bandwidth  dispersion 
issues,  let  eia°‘  and  f(t)eJ"J°‘  refer  to  the  transmit 
waveforms  in  the  narrowband  and  wideband  situation 
respectively  across  an  N  element  array  employing  M 
pulses.  Here  (see  Fig.  1) 


representing  the  spatial  steering  vector  along  Qk  and 
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Fig.  1  Narrowband  and  width  band  transmit  signals. 


for  the  l1',  2nd ,  ■■■  M'h  pulse,  we  get  the  space-time  vector 
z(t)  due  to  a  single  scatter  to  be  [2]  (details  omitted) 


zi(0 


z  (0  = 


(14) 


where 


Let  yf\t)  =  f(t)ejc0°‘  represent  the  scattered  return  at  the 

reference  sensor  due  to  the  first  pulse  from  a  ground 
location  at  an  azimuth  angle  0k  for  a  specific  range.  In  that 

case,  the  i'h  sensor  output  (suppressing  the  noise  component) 
is  given  by 

=  1  =  1,  2,  -  N  (8) 

where 

T  dsmfy  (9) 

1  c 


=  (15) 

This  gives 

z(t)  =  fk(t)os(0k,  ojdt)  (16) 

where  s( 6k ,  0Jd  )  is  as  defined  in  (2)  -  (5)  and  °  represents 

the  element-wise  (Schur-Hadamard)  multiplication.  Here 
fk  (/)  represents  an  MN  x  1  transmit  signal  dependent  vector 

whose  (iN+n)"‘  element  is  given  by  (see  (13)) 


represents  the  interelement  time  delay  corresponding  to 
azimuth  angle  Q  .  This  gives  the  first  pulse  output  across 

the  array  to  be 

yjt)  =  [f(t)eiwJ,  \)rx)eJ<,<N-^ .  (10) 
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Let 


Tk(T)  =  E{fk(t)rk(t  +  T)}>0  (18) 


Similarly,  the  received  signal  at  the  reference  sensor  due  to 
the  n,h  pulse  is  given  by 

y«\t)  =  f(t-{n-\)T2y<<-^  (11) 
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where  ® ,  is  as  defined  in  (5)  and 
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represents  the  interpulse  delay  for  azimuth  angle  0k .  After 

. th 

demodulation,  the  output  at  the  1  sensor  is  given  by 


zf(0  =  -1)t2 


-(i-  l)r1)e“M',“1K‘  , 

i  =  1,  2,  •  •  •  N,  n  =  1,  2,  •  •  •  M. 


(13) 


Notice  that  the  time-delayed  version  of  the  baseband 
waveform  in  the  return  signal  is  characteristic  of  the 
wideband  situation.  Using  (11)  -  (12)  with 

z  (t),  n  =1,  2,  M  representing  the  array  output  vector 


represent  the  covariance  matrix  associated  with  the  return 
wideband  signal  vector  ft(t)  in  (16).  We  will  assume  that 

the  actual  return  signal  components  in  (17)  resulting  from 
random  scatter  modulations  are  stochastic  in  nature  with 
wide  sense  stationaiy  behavior,  such  that  its  power  spectrum 
coincides  with  the  original  wideband  transmit  signal 
spectrum  in  Fig.  1  (b).  This  is  a  reasonable  assumption, 
since  the  returns  are  inherently  stochastic  due  to  the  physics 
of  the  problem,  nevertheless  their  spectral  content  is 
assumed  to  be  dictated  by  the  transmit  signal.  As  a  result, 
the  autocorrelation  function  TA  (r)  in  (18)  is  dictated  by  the 

inverse  Fourier  transform  of  the  baseband  version  of  the 
wideband  signal  power  spectrum. 

Extending  the  summation  over  all  scatterers  in  (16),  the 
array  output  covariance  matrix  in  the  fully  wideband  case 
takes  the  form 

Rz  =E{z(t)z\t)}  =  £/>*  Tk  os X  +a-\  (19) 

k 

where  T/;  =Tt.(0)  •  For  example,  a  low  pass  power  spectrum 
for  the  transmit  signal  gives 

Tk(m,  n)  =  sincB0((nl-n2)T2(k)  +  (i1-i2)Tl(k))  (20) 
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where  nt,  n2,  i]  and  i2  satisfy 


(26) 


m=ixN  +  n1,  n  =  i2N  +  n2.  (21) 

In  the  special  case,  such  as  when  the  azimuth  spread 
generates  only  a  small  angular  dispersion  A  9  due  to  a 
strong  main  beam  pattern,  the  Tt  in  (19)  may  be  replaced 

by  a  positive  definite  matrix  T  that  is  valid  for  all  k .  In 
that  case  (19)  reduces  to 

Rz=ToR^  (22) 

where  R  represents  the  narrowband  situation  in  (6).  It  is 

well  known  that  even  for  moderate  bandwidth  the  Schur- 
Hadamard  operation  in  (22)  results  in  spreading  the 
eigenvalue  spectrum  of  R  thereby  increasing  the  clutter 

degrees  of  freedom  [1],  This  is  shown  in  Fig.  2  for  a  14 
sensor  16  pulse  array  with  CNR  =  40  dB  and  T/;  as  in  (20). 
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Fig.  2  Clutter  eigen  spectra  for  narrowband  and  wideband  cases. 


III.  Wideband  analysis 

To  fully  understand  the  implication  of  the  Schur- 
Hadamard  product  operation  in  (16)  -  (19),  let  us  consider  a 
simplified  version  of  (19),  where  the  corresponding 
narrowband  scatter  set  corresponds  to  a  single  scatter  from 
90 ;  i.e„ 


Kx  =  P,MO„,ojdy(d0,ojd)  (23) 

and  hence  from  (19) 

R-=T0o  {P0s{90,cody{90,codo)).  (24) 

To  make  further  progress,  let 

MN 

T  =  A{  e,  e*  >  0  (25) 

1=1 

represent  the  eigen  decomposition  of  T  in  (18). 
Substituting  (25)  into  (24)  we  get 
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where 

p*  =  pA>  o 

and 
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(27) 


(28) 


Observe  that  (26)  represents  a  bunch  of  MN  uncorrelated 
returns,  all  of  them  associated  with  the  single  scatter  located 
along  0o  .  Interestingly,  s{  in  (26)  represents  a  spatio- 

temporally  amplitude  modulated  steering  vector  associated 
with  the  Doppler  frequency  a>d  and  location  Qg . 

Let  dA  represent  the  ordinary  DFT  vector  associated  with 
the  eigenvector  e;  •  The  entries  in  e,  correspond  to  a 
double  sampling  period  of  r,  followed  by  t2-  Asa  result,  it 
can  be  shown  that 


e 


k 


MN 


(29) 


where  dk(n)  represents  the  n'h  entry  of  d,  and  c0  is  a 
normalization  constant.  Here  i  and  j  represent  the 
unique  solution  to  the  equation  n  =  (jrl-l)N  +  in  , 
j„  =  1,  2 ,-M,  i„=  1,  2 ,-N. 

Using  (29)  in  (28)  we  get  the  modified  steering  vectors  in 
(26)  -  (28)  to  be 
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where 
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From  (30)  -  (31),  the  amplitude  modulated  steering  vector 
s  in  (26)-(28)  represents  a  coherent  sum  of  steering  vectors 

originating  from  locations  Q  n  =  1  — >  MN  with  Doppler 
components  a>d  ,  n  =1  MN  as  in  (31).  Referring  back  to 
(26),  for  each  scatter  return,  the  effect  of  finite  bandwidth  is 


3 


to  introduce  an  apparent  jittering  effect  both  on  angle  and 
Doppler  domains  so  that  a  bunch  of  uncorrelated  returns  are 
generated.  Each  such  uncorrelated  return  contains  a  set  of 
coherent  returns  with  different  directional  and  Doppler 
components  as  in  (31).  This  is  illustrated  in  Fig.  3. 
Observe  that  the  vector  sk  in  (30)  contains  several  coherent 

returns,  and  it  does  not  correspond  to  any  physical  steering 
vector.  As  a  result,  the  adaptive  processor  based  on  (26) 
will  not  be  able  to  null  out  either  the  original  return  from 
location  0  or  any  of  the  coherently  combined  vectors. 

From  (26),  the  wideband  data  z (t)  in  (16)  can  be 
equivalently  expressed  as 

MN  MN  MN 

*(t)  =  X  ak  Z  di  (”)  S(0r,  ,  )  =  Z  “A  (32) 

k= 1  n- 1  k-\ 


and  it  represents  MN  uncorrelated  bunches,  each  bunch 
containing  MN  coherent  returns  from  various  jittered 
locations  around  with  ()  ■ 

o 


Uncorrelated 
bunches  of 


(a)  Narrowband  case  (b)  Wideband  case 

Fig.  3  Narrowband  returns  and  effective  wideband  returns. 


Fig.  4  shows  the  SINR  loss  of  the  adaptive  processing 
with  and  without  bandwidth  dispersion.  From  there,  the 
clutter  null  is  wider  due  to  the  wideband  operation  nature  of 
problem. 


Presently,  there  are  only  two  basic  classes  of  techniques 
that  have  been  proposed  to  remedy  the  wideband  STAP 
clutter  cancellation  problem:  Subband  based  methods  or  3D 
STAP  [3]-[4] .  The  subbanding  approach,  as  the  name 
implies,  attempts  to  transform  the  wideband  processing 
problem  into  a  parallel  set  of  narrowband  2D  (angle- 
Doppler)  STAP  filters.  Since  a  single  2D  STAP  processing 
thread  is  already  a  quite  substantial  processing  burden,  a 
bank  of  say  L  subband  processors  would  represent  an 
L  x  100%  increase  in  processing  burden! 

The  situation  for  3D  STAP  is  even  worse.  “3D”  here 
refers  to  including  the  fast-time  (i.e.,  reciprocal  of  receiver 
bandwidth)  or  equivalently  instantaneous  frequency  domain. 
In  other  words,  the  dimensionality  of  the  STAP  problem  is 
increased  from  NM  to  NML,  where  N  is  the  number  of 
spatial  degrees-of- freedom  (DOFs),  M  is  the  number  of  so- 
called  “slow-time”  or  Doppler  DOFs  [1],  and  L  is  the 
number  of  fast-time  or  instantaneous  DOFs.  Since 
processing  burdens  are  roughly  on  the  order  of 
0(/fc2)  — M9(k3) ,  3D  STAP  can  result  in  an  L2  — >•  L3  -fold 
increase! 

The  above  processing  burden  increases  are  in  addition  to 
the  many  other  practical  difficulties  associated  with 
implementing  such  adaptive  filtering  schemes,  such  as 
sufficient  sample  support  for  weights  calculation  [1]  -  a 
problem  which  is  already  acute  for  2D  STAP. 

A  practical  framework  for  addressing  the  above  issues  is 
hierarchical  processing  in  which  an  “energy  preserving” 
lower  resolution  processing  stage,  to  which  2D  STAP  could 
be  applied,  is  followed  by  wideband  processing  only  if  a 
detection  was  obtained  in  the  first  stage.  Even  for  the  most 
stressing  MTI  applications,  detections  represent  a  very 
small  fraction  of  all  possible  resolution  cells — thus  resulting 
in  a  substantially  reduced  wideband  processing  burden. 
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